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. Continuous One-dimensional models supporting extended states are studied. These delocalized 

' states occur at well denned values of the energy and are consequences of simple statistical correlation 

^\ , rules. We explicitly study alloys of 5-barrier potentials as well as alloys and liquids of quantum wells. 

The divergence of the localization length is studied and a critical exponent 2/3 is found for the 3- 
^q' barrier case, whereas for the quantum wells we find an exponent of 2 or 2/3 depending on the well's 

parameters. These results support the idea that correlations between random scattering sequences 
break Anderson localization. We further calculate the conductance of disordered super lattices. At 
the peak transmission the relative fluctuations of the transmission coefficient are vanishing. 

O ' 

■ PACS numbers: 73.20.Jc, 73.20.Dx, 73.40.-c, 73.61.-r, 71.30.+h, 71.23.An, 71.23.-k 
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INTRODUCTION 



Recently the interest has increased substantially in ordcij-to find theoretical evidence for the breaking of Anderson 
localization by internal correlations in disordered systemsEni3. For instance, in ref.111 a tight-binding Hamiltonian with 
site correlations between the diagonal and the off-diagonal potentials was studied. A divergence of the localization 
length was obtained with a critical exponent v — 2 outside the band edges. In Refs-cTQ a simple discrete model 
' with correlations expressed by pairing (dimer) was studied. Their main result was the existence of a divergent 
localization length at some critical energies. In fact, this divergence is related to the existence of delocalized states 
found explicitly in ref.Q, which can be related to random phases and thus to extended states. Similar dimer models 
O | for classical systems were studied-,ip-jref.El. Using the inter-relation between a disordered Kronnig-Penney system and 
the dimer model (Poincare map)B~H3 an infinite set of such delocalized states exist. The generalization of the dimer 
y—( \ modeL-tp the N-mer case was considered, for instance, in refs£3ll3. In this context perturbative methods were used 
> ! in refill 

In this way, delocalized states have been found 'contradicting' the usual belief that all eigenstates are localized in 
one-dimension (1-d). In fact, this is only apparent because_the theorems related to localization in 1-d consider strictly 
uncorrelated random potentials Xaee for instance Ref-LETcHI). Interesting numerical simulations were carried out for 



binary disordered systems in ref.E3. 



■ Usually these correlated disordered systems support delocalized states for well defined energies. In this respect, a 
Q\ | 'band of delocalization' can be defined when the length L of the sample is smaller than the localization length L C (E), 
which diverges near the critical energy E c . Quantitatively, the region of delocalization can be defined approximately, 
as 

S3 ; 

where v has been determined explicitly in some models. Evidently this definition of a 'band of delocalization' is 
O ■ controversial because it depends on the existence of states in the interval \E — E c \ and, moreover, is rYSiy different 
, of the usual one related to true delocalized (unnormalizable) states. For instance, in the dimer modelaa a band of 
delocalized states is found to follow as AE ~ y/L, where L is proportional to the number of impurities TV, nevertheless, 
the total number of states grows as L and therefore the relative number of delocalized states tends to zero as 1/y/L. 
It has been furthermost recently shown that additional perturbations on the dimer potentials lead to a suppression of 
these extended statesEil. Therefore these delocalization properties become more relevant in finite systems, as presented 
in the last part of this article. In quasi-one-dimensional systems, like models with dimerized interchain couplingst3, 
leading to a cantor-set-like conducting band, the finite temperature transmission can be greatly enhanced by the 
resonance energies. 
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As noticed, internal correlations break strong localization. In this paper we consider some simple models showing 
delocalization properties. For a general point of view we consider the one-dimensional disordered Hamiltonian of a 
particle in a random potential 

H=^+J2v l (x-x l ), (2) 
2m * — ' 

I 

n f/0 -ai<Kai , . 

V ^ X >-\=0 otherwise ' {6) 

thus, a sequence of scattering barriers centered at the random position xi without overlap. Every barrier is symmetric 
with random support 2a/ and the distance between two barriers is given by 

Xi + i - xi = di, (4) 

which is a random quantity. 

Between consecutive barriers the particle propagates freely and its wavefunction can be expressed as 

i)l(x) = Aie lkx +Bie~ ikx for x t + ai < x < x l+1 - a l+1 , (5) 

where k is the wavenumber related to the energy by E — h 2 k 2 /2m. Theoretical group arguments&il relate the 
amplitudes (^4/+i, -Bz+i) to (Ai,Bi) by means of 



Ai+i \ ( e- lkx ' \ / m ftW e ikx ' \ / M 
B l+1 J I e lkx ' J { /3f a? J I er lkxi ) \B l 



(6) 



where the elements a, (3 of the collision matrix are inter-related like \a\ 2 — \(3\ 2 — 1 and we further on assume /3* = —j3 
(symmetric barrier). This form of the collision matrix is related to the fact that the time- independent Schrodinger 
equation is real. 

Defining the wavefunction just after the collision with the I— th barrier as 

tP+ = Aie lk(xi+ai) + S z e- lfc(:E!+Q!) , (7) 

a condition of derealization is 

^+1 = ±^ + , (8) 
which relates the elements of the collision matrix as 

a[e ik(d l +a l + 1 ~a l ) _|_ o* e -ik(dt+a l+ i+ai) _ ^ 

The above delocalization condition is a sufficient one and, as expected, is not verified for totally uncorrelated systems. 
Condition (9) yields using \a\ 2 — \/3\ 2 — 1 the following simple delocalization condition 

H = l and 01 = 0, (10) 

which leads to \Ai +1 \ = \ A t \ and \Bi +1 \ = \B t \. 

Condition (9) inter-relates the random parameters like a/ and di in the disordered system and is therefore a source 
of correlation between these random parameters. In fact, in some simple systems, relations (9) and (10) can be 
satisfied and we present two examples: The <5-barriers-alloy and the quantum well liquid-alloy with interrelated 
random parameters. 

To end this section, we remark that any relation similar to (6) can be formally written as a tight binding Schrodinger 
equation (Poincare map) by considering 

ai = aie lk{dl+ai+1 ~ ai \ fa = ^ e ^(rfi+«i+i+«0 (H) 

and 

+ A^ti = V^+, (12) 

where 

D(k) 

Dl = — / ' and V5 = {A+i(fii+#)+A(&f-i-#-i)}- (13) 

liylmoii — lmpi) 

We notice that the above expression is only a formal one. For models such as the ^-barriers sequence, Di is well 
defined for any k with the arbitrary choice D(k) — 2isinfc. Nevertheless, in other systems singularities for some 
values of the wavenumber k can appear. 
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I. DELOCALIZATION IN A 5-BARRIER SEQUENCE 



Consider a sequence of 5-potentials which are statistically distributed over lattice sites (alloy). Thus consider the 
sequence 

Vi(x) = Vi5(x-xi) and d l = d, (14) 

where Vi are random uncorrelated parameters and d is a constant lattice parameter. In this case, the elements of the 
collision matrix are given by 

Vi Vt 

a; = l + i— and (3 l =i—. (15) 
Ik Ik 

The derealization condition (9) becomes explicitly 

e^ferf _ ^ sin(fcd) = ±1. (16) 

In general, for any arbitrary momentum k, this condition does not hold because of the random quantity Vi. Neverthe- 
less, if k = mr/d ( n 6 Z*), as first observed by IshiiE£l, we have a set of dclocalized states where \ipi+i\ — \tpi\- This is 
at first sight very surprising because of the simplicity of the disordered model. Evidently, correlations exist and are 
related to the choice di — d for any In fact, assuming a sample with TV barriers, the usual 27V random parameters 
characterizing the uncorrelated system are reduced to N because of the constraint di = d. 

In this model the physical interpretation is simple, in fact as long as k = mr/d the electron does not 'feel' the 
random potential because the distance d, between consecutive barriers, is a multiple of its wavelength. 

At this point, it is interesting to study the divergence of the localization length L C (E) near to the critical energy 
= 2m(^T") 2 - Using the Poincare map (12) for this model we have 

Vi 

ipi+i + ipi-i = 2{cos kd— sin kd}tpi. (17) 
For e = E c — E <C 1 and H/2m taken as unity, equation (17) can be rewritten 

+ Vi-i = ± ^2 - Vje^jW ( 18 ) 

This last model was extensively studied in the limit e <C 1 by Derrida and GardnerEl. They calculated the complex 
Lyapounov exponent 7, where the real part corresponds to the inverse localization length and the imaginary part to 
7r times the integrated density of states. Their results can be expressed as follows: 



7m( 7 )~X 2 e 2 /3( V; 2)i/3 



(n7r) 



2/3 ; 



where K\ = 0.29 . . . and K2 = 0.16 . . . and (•) is the average over all impurities. From (19) it is straightforward that 
the inverse localization length L" 1 scales as 

L- 1 ~e 2 / 3 (V 2 ) 1/3 (20) 

and the density of states is 

p(e) = djm-f(e) - e" 1 / 3 . (21) 

The main result of this section is to show that a continuous disordered model where the potential amplitudes are 
random but placed on a lattice, exhibits extended states for discrete values of the energy. We recall that for this 
model the result was first pointed out by Ishiilij. However, the exponent 2/3 we obtain is different to the 1/2 obtained 
by Ishii, who considered a Cauchy distribution. The correlation in this model comes from the fact that the impurities 
are placed on a regular lattice. This in fact is a very natural physical assumption, as one can imagine an alloy where 
different atoms are placed randomly, described by the random amplitude of the potential, but on a lattice. This 
system would lead to a similar model studied above. In fact the models where the locations of the impurities are 
continuously random, no extended states have been found so far for delta like potentials. For finite potentials they 
do exist as we will mention in the next section. 
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II. DELOCALIZATION IN A QUANTUM WELL SEQUENCE 



The case of delta impurities considered in the previous section is very special and cannot be used in order to 
describe more general impurity potentials. This is why it is of interest to study other types of potentials, which 
could eventually describe more general impurity potentials. One of the simplest examples of such a potential is the 
rectangular well potential. In principle one could approximate any potential with rectangular wells if the elementary 
width is small enough. Therefore we will restrict our considerations to rectangular well potentials. 

As an important result we observe the existence of extended states for systems where we have continuous spatial 
disorder and discrete shape disorder. This has to be contrasted with the case of delta impurities where extended 
states exist only when we have shape disorder alone, i.e., where the amplitudes of the impurity potentials are random. 

The aspect and definitions of the rectangular well potential are shown in fig. 1: 
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FIG. 1. Rectangular well potentials 



For this case we have 



and 



ai = (2pik cos(2a lPl ) + i(pf + k 2 ) sm(2aipi)) (22) 



_ A {pj -k 2 )sm(2a m ) 

01 ~ l 2p7k ' (23) 



where pi = \J\k 2 — V\\, We have to point out that these results are only valid when E > Vi- For the case E < Vi we 
have to replace cos(-) and sin(-) in (22,23) by their hyperbolic counterparts. 

Using the derealization condition (10) we have an extended state when (3i — 0, thus 

? 9 

E = k 2 = Vi + -t-j- with ni an integer. (24) 

This means that if we have a discrete set of impurity configurations, we have critical energies E c at which the states 
are delocalized. The impurities have just to be constructed with the help of equation (24) once we replace E by E c , 
where we can take different values for Vi or ai by changing n; . A very remarkable aspect of these disordered models is 
that they are independent of the space disorder expressed by di as was already appreciated by TongtS for a different 
model. 

We evaluate the exponent for the following case: suppose di = d and a; = a, i.e., we are left with a pure shape 
disorder. Let k 2 = E c be the critical energy defining the correlation (24) of the disorder and writing the Poincare 
map (12) around k — k c + e, with e very small, we obtain 

ipl+l + in-i = {2 cos (k c (d - 2a)) + eft}^, (25) 

where is random and function of Vi- For the case where k c (d — 2a) /it is an integer we obtain the same result as 
for the delta-barrier case (18-21), thus v = 2/3. When k c (d — 2a)j-K is not an integer we can also use Derrida and 
Gardner'scj result for the case inside the band, yielding v = 2. The most general case, where di and ai are random, 
and when no particular configuration is used, we obtain numerically the exponent v = 2. 
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The interesting aspect of these continuous disordered systems is that they can describe a real physical situation of 
electrons in a certain potential. The problem with continuous models is that the expressions of the transfer matrices 
can be much more complicated except if one considers delta impurities or rectangular well potentials. For more general 
impurity potentials one can either approximate them with a sum of rectangular barriers or use the semiclassical WKB 
approximation. This last approximation can, however, have dramatic effects on the localization properties, especially 
around some critical energies, as they are very sensitive to small perturbations. 

Summerizing the localization length divergences for the (5-impurity and quantum wells we have: 

r ,_v 1 , J = 2/3 ^-barrier and k c (d— 2a) /ir an integer for the quantum well ,„_.. 

L ^ E > ~ \e-EJ" ' Where V = = 2 otherwise (26) 



III. FLUCTUATIONS AND CONDUCTANCE 



In one dimensional disordered systems, the relative fluctuations of the transmission coefficients diverge with the 
size of the systernES. Therefore the average transmission is a bad statistical quantity. However at the critical energies 
discussed above, the states are deterministic and therefore one would expect no fluctuations. In the following we 
evaluate the relative fluctuations explicitly and calculate the conductance for superlattices around these critical 
energies. 

Around the critical energy E c we evaluated the localization length dependence on energies, therefore T/v ~ e~ L ' L ", 
where T/v is the transmission coefficient after N barriers and N is proportional to L. The relative fluctuations are 
then given by 

Using the central limit theorem, valid for the logarithm of the product of random determinant one matrices, the 
relative fluctuations of the localization length can be expressed as 

^ ~ TV" 1 / 2 , (28) 
L c 

hence using (26-28) 

AT 

— ~NV 2 \E-E C \ V . (29) 

In usual one-dimensional disordered systems the average distance between eigenvalues is of the order JV _1 , therefore 
using \E — E c \ ~ N , we obtain 

AT 

— ~\E-E C \»-^. (30) 

This demonstrates that for the case of interest here, both cases v = 2 and v = 2/3 have vanishing relative fluctuations 
near the critical energies. This fluctuation analysis gives also a bound on the possible divergence exponent v. In fact 
as long as v > 1/2, the average transmission around the critical energies is a good statistical quantity. 

Many physical systems can be very well approximated with either rectangular well potentials or delta potentials. 
This is the case for superlattices of heterostructures. One can grow very precise, up to the atomic precision, layers 
of GaAs. Then one dopes differently one layer with Al and then again the same layer and one obtains in this way 
a single so called delta doped layer. This can then be repeated as often as one wishes. This last structure is very 
well described in the direction of the layers by delta impurity models. Instead of doping differently only one layer 
we can also dope a finite number of layers differently. In this case the potential looks more like a rectangular well 
potential, where the amplitude depends on the concentration of the doping used. In this way it is possible to grow 
one-dimensional quantum wires where both directions of the layers are etched, two-dimensional systems where only 
one of the direction of the layers is etched and finally three-dimensional systems. r— ■ 

The two probe conductance is simply given by the sum of the transmission coefficientsEil, thus 

G( M ) = J2 T ^ ■ X' (31) 

k\E=fi 
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where /x is the chemical potential. 

For a one-dimensional quantum wire the conductance is just given by the transmission coefficient, as there is only 
one conducting channel. Around the critical energies we evaluated the localization length dependence on energies, 
therefore using TV ~ e~ N / Lc , where T/v is the transmission coefficient after TV barriers, we obtain conductance peaks 
at the critical energies of the form 

G(n) ~ e-rlM-Enl^ (32) 

where fi is the chemical potential, 7 depends on the size and disorder of the system and E n are the critical energies. 
For the delta doping E n = -^(^j 1 ) 2 and v = 2/3. For the rectangular case we have one critical energy E c and the 
random Al doping potentials are given by 

n 2 TT 2 

Vi = E c with m a random integer (33) 

and a is the width of the doped layers. The width di of the undoped layers can be random and we obtain v = 2. For 
the special case where di = d is constant the value of the exponent v depends on the ratio E c (d — 2a) 2 /it 2 . If this 
ratio is an integer the exponent is v = 2/3 but remains v = 2 otherwise. The interesting feature of this system, is 
that it would behave as a perfect filter. The band-width could be controlled by the size or 7 and the desired critical 
energy by the choice of dopants. 

In order to evaluate the conductance of the three-dimensional system, as defined above, we have to sum over all 
possible channels. In the directions (x,y) perpendicular to the grown layers the density of states is simply the usual 
box filling density. In the infinite limit we have to integrate over the Brillouin zone which yields, 

G~ / e -i\n-kl-k 2 y -Ec\> dkxdky (34) 

The plots of the conductance from eq. (34) as a function of the chemical potential are given in fig. 2 for different 
system sizes characterized by 7 and for the two cases v — 2 and v — 2/3. 

The two-dimensional case is presented in fig. 3. In this case the behavior is consequently different for v = 2 and 
v = 2/3. Indeed for v = 2/3 we have a singular derivative, which could be observed experimentally, in order to 
distinguish between the two cases. 




Chemical potential Chemical potential 



FIG. 2. Conductance for v = 2 and v — 2/3 (right graph) as a function of the chemical potential, for the case where the 
resonance energy is E c = 1. The different curves represent different values of 7 = 2 n , with n = 1,2,3,4,5. The uppermost 
curve represents the case 7 = 2. 



6 





Chemical potential 
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Chemical potential 



FIG. 3. Conductance for the two-dimensional case with the same parameters as in fig. 2. 

The quantitative values can also be estimated. In the delta doping case the typical Fermi wavelength is of the order 
of 200A and the typical layer thickness is of the order of 5A. This means that we have to construct delta doped layers 
with around 40 layers in between. In order to obtain the disorder we just have to chose randomly the dopant or the 
concentration of the doping for the delta layers. The same reasoning can be held for rectangular potentials, where 
one just has to chose a particular potential configuration. 

CONCLUSIONS 

In the first section we demonstrated the existence of extended states in continuous disordered systems with delta 
impurities placed on a lattice and obtained for the exponent describing the divergence of the localization length 
v = 2/3. This model hence demonstrates the existence of extended states in systems with only shape disorder. In 
the second section, using rectangular potentials we demonstrated that models with shape and space disorder can also 
exhibit extended states. The only constraint on the disorder distribution is that it must be discrete for the shape 
disorder but can be continuous for the space disorder. In this case the correlation in disorder leading to the existence 
of delocalized states is simply the fact that the distribution mentioned above is discrete and obeys relation (24). 
Therefore one has to be very careful when one discusses disordered systems using discrete disorder distributions as 
there can exist singularities in the spectrum. The exponents of the localization length divergence are very interesting. 
In fact, depending on the parameters of the quantum wells we obtain different exponents, either v = 2/3 or v = 2. 

The relative fluctuations of the conductance around these critical energies are vanishing. Therefore the conductance 
is a well defined statistical average and represents the transport properties of these disordered systems. 
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